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This two volume report documents the v;ork accomplished during 
amended Phase II of Contract Number NAS 9-14636. The system 
tradeoff analysis and system performance study that led to 
three antenna pointing and acquisition algorithms for hardware 
simulation is provided in Volume IV. The construction of 
results obtained from the hardware simulation are described 
in Volume V. A high level software description and a detailed 
software documentation of the computer programs is to be 
found in Volume V. The first three volumes of this report 
documents the work accomplished during Phase I of this contract. 
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OBJECTIVES 

LINCOMs overall objective under this contract was to develop 
a hardware simulation which could emulate the Shuttle's Ku-Band 
Antenna Pointing and Signal Acquisition System. The desire was 
to develop a simulation in which the antenna pattern, the TORS 
search volume, the a priori probability distribution of TORS 
satellite position relative to the Shuttle, and the antenna scan 
procedure could be selected for the purpose of pre icting 
performance for various Shuttle/TORS antenna pointing and 
acquisition scenarios. 

The simulation was developed under the constaints of assuming 
a fixed network operation, the system must be real world imple- 
mentable, it must be cost effective, the program execution time 
must not be excessive, the modulation technique is PN/BI-41/BPSK, 
and optimum performance is desirable. 
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This document presents the results pertaining to the trade- 
off analysis and performance of the Ku-Band Shuttle antenna 
pointing and signal acquisition system. The study was performed 
assuming the existance of various antenna scanning trajectories 
and various signal acquisition algorithms. The square, hexagonal 
and spiral trajectories were investigated assuming the TDRS 
postulated uncertainty region and a flexible statistical model 
for the location of the TDRS within the uncertainty volume. 

The scanning trajectories, Shuttle/TDRS signal parameters and 
dynamics and three signal acquisition algorithms were integrated 
into a hardware simulation discussed herein and documented in 
detail in Volume V. The hardware simulation is quite flexible 
in that it allows one to evaluate signal acquisition performance 
for an arbitrary (programmable) antenna pattern, a large range 
of C/Nq's, various TDRS/Shuttle a priori uncertainty distributions 
and three distinct signal search algorithms. 

Based upon the data made available during this contract 
period, certain Ku-Band forward link signal threshold characteristics 
were studied. The antenna pointing and acquisition threshold is found 
to be less than 60 dB-Hz with an acquisition time dependent on the 
antenna scan procedure and acquisition algorithm implemented. 

Various techniques are discussed in this report and a computer 
program is presented from which these can be evaluated. It 
appears that there will be no problem in meeting the system 
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spec even in light of the uncertainty associated with the 

possibilities of antenna sidelobe acquisition. The details are 

provided herein. Assuming a single channel monopulse system and a 

30 MHz IF bandwidth, a monopulse tracking loop bandwidth of one Hz, 

the monopulse tracking jitter is 0.11 degrees at C/Nq = 54 dB-Hz. 

The Costas loop arm filter bandwidths can be chosen to be ten 

times the data rate such that no false-lock problem occurs 

during carrier acquisition. Assuming a loop bandwidth of 

= 3 kHz, the loop jitter is 10 degrees at C/N^ = 60.4 dB-Hz. 

Acquisition can be accomplished in less than 10 seconds. 

For the purpose of antenna scanning analysis, a reference 

coordinate system, whose z-axis is in line with the center 

axis of a specified uncertainty cone of the TORS position, is 

chosen for the relative relation between the TORS and the Shuttle 

antenna. The scan path of the Shuttle antenna can be projected 

onto the (x,y) plane of the coordinate system; while the uncertainty 

cone of the TORS position can be described by a circle in the 

(x,y) plane. The uncertainty in the position of the TORS is 

modeled by a truncated Gaussian probability density p(x,y) with 

2 2 2 

uncertainty parameters '7^ = = a . By changing the value 

of this variance parameter, the model is sufficiently general 
to include a uniform distribution of TORS position uncertainty 
(op=“) to one which specifies the position with probability one 

(ap=0). 

In the study of antenna scanning three types of trajectories 
are proposed ard evaluated in terms of average scan time and the 
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structure of coverage over the uncertainty cone of the TORS 
position. The analysis technique used for finding average 
scan time is discussed and an illustration of the technique 
for a specific type of scanning trajectories and motion of the 
Shuttle antenna is given. The results show that the spiral tra- 
jectories is in general better than the other two trajectories, 
especially for a constant velocity along a trajectory, since the 
path of a spiral trajectory (from the center of uncertainty cone 
to its edge) can be shorter than the other two. 

From the Ku-band system specifications, it has been 
recognized that the variation of the received signal level at 
the Shuttle (due to TDRS EIRP path loss and antenna pointing loss 
variations) varies as much as 23 dB. If the antenna sidelobes are 
not sufficiently suppressed, a potential problem called sidelobe 
acquisition may cause the degradation in system performance. 
Therefore in this tradeoff analysis study, three different 
aroiiisition strategies are proposed and evaluated. The first 
acquisition strategy may be used for the case where the sidelobe 
acquisition does not impose a problem to the system. The other 
two strategies are primarily designed to avoid the sidelobe 
acquisition (especially the first sidelobe). 

In this study, we assume that the first sidelobe of the 
Shuttle antenna is 17.5 dB suppression from the main lobe of the 
antenna so that we can study the effectiveness of the acquisition 
algorithm proposed. Surely, one should note that the sidelobe 
acquisition can be overcome by tapering down the sidelobes of 
the Shuttle antenna below 23 dB. However, it has been 
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found that with one particular acquisition algorithm studied, 
a strong sidelobe of the antenna can speed spatial acquisition 
significantly. 

The scan schemes and acquisition algorithms are all integrated 
by a computer simulation program. This software package, discussed 
in Vol. V, provides a useful tool for predicting the performance 
of the Ku-band antenna pointing system. Various options are 
available for users to do tradeoff studies on system parameters 
in designing an acquisition system. The software package has 
been tested on UNIVAC 1100 Series computers and verified with 
analytical results. 
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I. THE TORS POSITION AND SHUTTLE ANTENNAS 

1 . Reference Coordinates 

In establishing the Ku-band communication link between a 
Space Shuttle receiver and the TORS, the Shuttle receiver must 
perform spatial search for the TORS position. Initially, it 
is assumed that the transmitter antenna of the TDRS illuminates 
the Shuttle with a spread spectrum signal and also that thr 
position of the TDRS is nearly stationary during the process 
of spatial search of the Shuttle antenna, since the TDRS 
dynamically is moving along a figure 8 trajectory slowly. 

To discuss the spatial search of the Shuttle acquisition 
system, we shall establish a reference of coordinate systems. 

In this study it is assumed that the size and center of the 
uncertainty cone within which one can locate *he TDRS, are known 
at each Shuttle receiver. The geometric diagram of Fig. 1 
depicts in general the relative relation between the Shuttle 
receiver and the TDRS. However, for the convenience of 
references we adopt a coordinate system whose z-axis is in line 
with the direction of the center of the uncertainty zone. Denote Az 
to be azimuth coordinate and El to be elevation. We shall project 
Az and El onto a plane perpendicular to the z-axis, as shown in 
Fig. 2, so that the uncertainty zone can be confined to a circle 
whose radial coordinate represents Elevation and angle 
coordinate Azimuth. Throughout this report, we shall use this 
projected plane to discuss the spatial search of the Shuttle 
antenna. 
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The position of the TORS relative to the Shuttle receiver 
is only known up to a cone of angle 9^. Nothing has been 
specified about the statistical distribution of the TORS over 
the cone. In order to characterize its distribution, a 
"truncated" Gaussian function is assigned to (x-y) coordinates 
of the uncertainty zone on the projected plane. The "truncated" 
here means no tail for the Gaussian function, i.e., the probability 
of finding the TORS at one point in x (or y) direction is of 
Gaussian form within the uncertainty region and zero outside 
the uncertainty zone. The value of variance Op of the Gaussian 
distribution reflects the statistical properties of finding 
the TORS around the neighborhood of the center of uncertainty 
cone. The smaller the value of Op, the higher the chance 
of finding the TORS near the center of the cone. For a higher 
value of Op, it approaches to a uniform distribution. 

2. Statistical Model of the TORS Positio n 

Now the probability distribution is expressed in terms of 
(x-y) coordinates in the projected plane of Fig. lb. The 
Gaussian distribution function is 


♦(x,y) = 


1 


2ira_ 


exp 


2 ^ 2 
X +y 


2o 


(1) 


2 


wher’ Op is the variance. The probability distribution of the 
TORS at (x,y) is 


f(x,y) = 


k$(x.y) 


2 ^ 2 2 
elsewhere 


( 2 ) 


where k is a normalized constant. One can show that 
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k = 


\ - exp 


(-4) 


by imposing the condition 
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'’age is 


f(x,y)dxdy = 1 . 


Note that the function f(x,y) has the same zero mean as that of 

2 

♦(x,y), but it has different variance from Op. One can find that 
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In practice, it is more interesting to know the statistical 
characteristics along the direction of elevation and azimuth. 
One can express f(x,y) in terms of polar coordinate system 
in the projected plane as follows 


/ kR 


f(R.e) = < 


2uo 


exp - 


0 < R < 0. 
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(5) 


el sewhere 


where R is for elevation and e for azimuth. It is obvious that 
the probability distribution in the azimuth is uniform over (0,2tt). 
Thus 



kfi exof.-i 

2 '■'M 2,2 
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elsewhere 


The mean and variance in the elevation can be shown to be 
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E(R^) = ko^p, 


!'-('* J-) (- 4)) 


9 2 

Variu.ice = E(R“) - mp 


(7) 


3 . Mean, Variance and Mo m ents f or the TDRS Position 

In order to show the relation between mean and variance 

2 

of the elevation distribution and the variance Op, a tabulation 

2 

of means and variances for different values of o is given in 

P 

Fig. 3. In addition, because of the need in analyzing the 
average scan time, the moments of elevation and azimuth 
distributions are also tabulated in Fig. 4 (all normalized 

4 . Shutt l e Receiver Antenna Pattern Model 

It is assumed here that the pattern of the Shuttle antenna 
has only one single pencil beam with the first sidelobe down 17.5 
dB from the main lobe. The beam width shall be 

referred to as the 3 dB contour of the antenna pattern. The 
type of the antenna pattern assumed here can be characterized 
as follows: 

J^(2a cos 0 ) 


G(o) 


a sin 6 


( 8 ) 


where 


0 = angle offset of the boresight axis of the Shuttle 
antenna from the direction of the TDRS position 


a = a parameter determined by 3 dB beamwidth of the 
antenna Oc 
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Figure 4. MEAN AND MOMENTS OF ELEVATION AND AZIMUTH 
DISTRIBUTION OF TDRS POSITION. 


Elevation: 

(I-l)TH CENTRAL MOMENTS 

1 l.OOOOOOOE+00 

2 3.1307825E-01 

3 2.6646434E-02 

4 2. 6520391 E-03 

5 2.2308420E-03 

6 6.2739134E-04 

7 3.5791496E-04 

8 1.5581140E-04 

9 8.471 5485E-05 

10 4.4132538E-05 

11 2.4764277E-05 

12 1.40-0367E-05 

13 8.1651569E-06 

14 4.8283923E-06 

15 2.9042068E-06 


Azimuth: 


MEAN, AND 
MEAN= 5. 


10 0 . 
11 8 . 
12 0 . 

13 1. 

14 0. 

15 4. 


CENTRAL MOMENTS 

OOOOOOOE-01 

3333333E-02 

2500000E-02 

2321429E-03 

3402778E-04 


8778409E-05 
8780048E-05 
06901 04E-06 
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J^(x) = Bessel function of first-order 

Note that here we assume the antenna pattern is spherically 
synmetric with respect to the boresight of the antenna. The 
function 2lJ^(x)l/x is plotted in Fig. 5. From the function one 
can see that the 3 dB point of the pattern is located at x = 1.61 
and that the first sidelobe peak is 17.5 dB down from the main 
peak and occurs at x = 5.1. With this information one can 
design the antenna with 3 dB beamwidth by choosing the 
parameter a to satisfy the condition 


where 


. 1.61 0.805 


0 g = half beamwidth (in rad.) 


(9) 


In what follows we set a-57.65 and this gives rise to a beamwidth of 1.6®. 
The antenna pattern used in this study is plotted in Fig. 6 with an amplifi- 
tion factor of 3 in angle and in Fig. 7 in rectangular form. Note that 
the first null point of the main lobe is at = 1.78° and that 
the angle of the first sidelobe peak is 0^ = 2.5°. 

The function G(e) gives the loss in the received signal 
power when the angle 0 is nonzero. The e, an offset angle, is an 
angle of the boresight axis of the Shuttle antenna off from the TORS. 

To compute 6, one should know the relative geometric relation 
between the TORS and Shuttle receivers. Using spherical coordinate 
system, we denote A to be a fector pointing toward the TORS 
with unit magnitude, say (1, B be a unit vector 

in line with the boresight axis of the Shuttle antenna, say 
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Figure 3. The Function 
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(1 ,Saz»S E|_)- Then one can express the offset angle 6 to be 


cos 0 = 


■ aTTb’ 


where 


lx| = length of a vector X 
• = scalar product of two vectors 


In scalar form the above equation becomes 


cos 0 = cos(aAz-SAz) “el ®EL ^ “EL ^ E’ 


If the TORS happens to be at the center of uncertainty region, 
i.e., a^z = 0 and = 0, one immediately has 

cos 0 = cos 6 el ® ^ EL 


Eq. (10) will be used in general to compute the offset angle 0 
for a given position of the TORS and the direction of the 
boresight axis of Shuttle antennas. The latter will be 
represented as a point in the projected plane as indicated in 




II. SPATIAL SCANNING OF SlJ.'TTl r AN’i ENNA 


1 . Scanning Trajnctoric .s 

In order to search spatially for the TORS position systematically, 
one should assign a regular path for the Shuttle antenna to sweep 
across the uncertainty region. Projecting ihe path, along 
which the Shuttle antenna sweeps through, onto the plane perpendicular 
to the 7-axis, we have a trajectory for the tip of boresight of 
the antenna and footprints for the 3 dB contour of the antenna 
pattern. To guarantee a complete coverage of the uncertainty 
region with a 3 dB contour footprint, one should properly design 
the trajectory. In this study, three kinds of scan trajectories 
are studied, viz. , 

(1) Square Trajectory. 

(2) Hexagonal Trajectory. 

(3) Achimede Spiral Trajectory. 

The characteristics of each trajectory will be discussed separ- 
ately. Here we'll point out that the first two trajectories 
will have uniform overlaps between any two 3 dB contour of the 
antenna pattern. The whole uncertainty area is divided into 
cells inside which every point can be illuminated by at least a 
3 dB contour of the Shuttle antenna. On the other hand, the 
third trajectory will have nonuniform overlaps, especially 
those cells near the center of the uncertainty zone. The 
overlap used here is defined to be the ratio of the portion 
o' the overlap of two cells whose centers on a trajectory are 
aligned in the same radical direction, to the 3 dB beamwidth 
of the antenna, as shown in Fig. 8. 
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20r 


( 11 ) 


Here we also assume that all trajectories start at the center 
of the uncertainty cone. 

(a) Square Trajectory 

As the name implies, the trajectory is square-shat- ed and 
the uncertainty area is divided into square cells with diagonals 
20g, as shown in Fig. 9. The center of each cell cain be easily 
determined from the geometry of the trajectory. One can see that 
the disadvantage of the square trajectory is chat some of the 
cells are outside the specified uncertainty cone. Obviously 
one can easily determine the number of terms of the trajectory 
and the number of colls needed in order to completely cover the 
whole uncertainty cone. 

Let be the angle of half i.he uncertainty cone, 2d the 
separation between two square cells, 2d = v^ 0 g and N the number 
of turns. Then from observation of geometry, one can easily 
express 


0„ + 2Nd > 0„ 
o — m 


or 


N > 


m b 


- 7? (^- ') 

This indicates that it takes about ^ l^urns of 

square trajectory to sweep through the whole area. One may 
also estimate the number of cells needed 


( 12 ) 


= 1 + 8 + 2x8 4 ...+ 8xN 
c 

= 1 4- 4N(N4-i) 


(13) 
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For instance, = 0.8'’, = 10*', then one needs 


N > 8.13 at least N = 9 turns 


and the number of cells 


N = 1+ 4 X 9 X 10 = 361 cells 
c 


Finally, one should note that the overlap for the square trajectory 
is 29.3%. 

( b ) Hexagonal Tra jectory 

Hexagonal trajectory requires less overlap between cells. 

The intersections with all neighborhood cells are the six vertices 
of hexagons which form a honey comb. The separation 2d between 
two hexagonals is *30g. as shown in Fig. 10. Obviously, the 
centers of the cells can also be determined easily. Here we 
can also estimate the number of turns required to sweep the 
uncertainty area to be as 




one can solve for N 


N 


ORIGINAL PAGE IS 
OP POOR QUAUTY 


m R 
2d~ 


1 



- 1 


(14) 


T7 


The number of cells can also be computed 


or 


= 1+ 6+2x6+. .. + Nx6 

c 


= 1 + 3N(N+1) 


(15) 


Comparing eq. (13) with (15), one can see that the number of 
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cells required for the hexagonal trajectory is less than those 
for the square trajectory, say Oo - 0.8'’ and o„ = 10'’* the 

D Ml 

trajectory needs 

N 7.1 at least N = 8 turns 

and 

N, = 216 cells 
c 

which is snaller than that for the square trajectory. 

(c) Archimede Spiral Trajectory 

When the boresight axis of the Shuttle antenna sweeps along 
a trajectory, the foot prints of a 3 dB contours of the antenna 
cover a uniform strip of area along the path. In order to 
utilize this fact, the Archimede spiral is most suitable because 
of its equal separation between trajectories. The Archimede 
spiral can be characte* ized as 

El = ” Az , or = koi 2 

where 

2d = separation between trajectory 

Note that the variable 02 (or Az) in the expression should not 
be limited to a range of 2ir. According to the definition of 
overlap eq. (11), the seraration d should be determined once 
the overlap of cells is specified. 

Unlike other two trajectories, the spiral trajectory lacks 
geometrical symmetricity. So how to choose the center of cells 
to ensure a complete coverage of the scanned area becomes a 
nontrivial problem. 

To explain the scheme of finding centers of cells along the 
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spiral, we shall start the first rell of the scan at El = 0 and 
Az = 0 and proceed along the spiral. We should Impose a 
condition that the circle of the second cell passes through 
the Intersection of the circle of the first cell and an 
auxiliary spiral. This auxiliary spiral falls In the Archimede 
spiral trajectory and can be characterized as below 


0. 


a- +d 

n c 




as shown In Fig. 8. 

Let A be the intersection point of the first cell and 
the auxiliary spiral, ® center of the 

second cell (El 2 .Az 2 )> Using eqs. (16) and (17), one can solve 
for Eljjj^, Az^^ and EL 2 .A 22 In sequence. 


El„, = Ojsnd 




and 


Az 

m 


1 ■ > 


El, 


- Az, + d 

It 2 


(18) 

(19a) 

(19b) 


The relation (19a) Is derived from the fact that the distance from 
point A to B Is 0g. Combining two equations (19), one can 
Iteratively solve for Az 2 » then El 2 * Repeating this procedure, 
one can find the centers of all cells to cover completely the 
scanned area. Here we shall formulate the general method. 

Let the center of the k— cell be known (Elj^.AZj^). To 

A ^ 

find the center of the next cell ((k+1) -) one should find the 
intersection of the circle of the k~ cell and the auxiliary 
spiral curve first. Denote the intersection to be (^^nik*^^mk^‘ 
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Then one can have the relation 


Az 


mk 


= AZj^ + cos 


■)/. \ 


( 20 ) 


to solve for Using eq. (17) one can compute After 

(k+l)» 


obtaining the values of Az^j^j^ and one can solve for Az 


iteratively, by the relation 


Az 


(k+1) 






®B ~^^ ^(k-H) 


2E1 


(k+1) 


llW 

E'mk ) 
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Similarly, one can obtain by eq. (17). 

From the above discussion of finding centers of all cells 
one can realize that it is rather complicated to analytically 
evaluate the number of cells required or number of turns needed 
to cover whole area of uncertainty cone. However, a computer 
program implementing the above scheme has been developed and can 
provide those answers for a specified uncertainty cone. 

(2) Motion of the Shuttle Receiver Antenna 

Once a trajectory of spatial search is chosen for Shuttle 
receiver antennas, one needs to determine what type of motion 
should be used in order to meet aspects of the total search time. 
Here, it is assumed that the receiver antenna can be moved 
either in a constant angular velocity or constant speed in any 
trajectory or a combination of them. The constant angular 
velocity can be characterized as 

do2 

= constant 


( 22 ) 
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and the constant speed along a trajectory is defined to be 


gr = constant (23) 

la long trajectory 

where s is the arc length along a trajectory. 

How should we choose the type of motion for the receiver antenna 

depends on the flexibility of gimbol systems, which drive the 
antenna to sweep, on the adoption of scanning trajectories and 
on the consideration of dwell time for each search cell. One 
important parameter in the design of acquisition system is the 
total scan time. Here we shall discuss the effect of types of 
motion on the scan time. We shall neglect the time required 
for signal detection since it is normally much shorter than 
the time for mechanical motion. 

Suppose that the beamwidth of the Shuttle antenna is 1.6®. 
the scan times to sweep over an uncertainty cone of 20® are 
plotted in Fig. 11 for constant angular velocity of Shuttle 
antennas with the spiral trajectory. For instance, 'or a constant 
rotation speed of 100“/sec, it will take 32 seconds to sweep the 
whole uncertainty area if the overlap of the spiral trajectory 
is 30X. To meet the required scan time, 55 sec, it is necessary 
for the Shuttle antenna to rotate at least 50®/sec, for 15% of 
overlap spiral trajectory. This is equivalent to 0.131 cycle/ 
sec. The scan times for other two types of trajectories are 
also shown in Fig. 11. 

On the other hand, the curves given in Fig. 12 show the 
scan time for constant speed along the trajectory. Note 
that the unit of speed is the same as the unit in elevation 


I I 






per sec. The uncertainty cone is also assumed to be 20®. The curves 
show that the Shuttle antenna requires tc move 4.9°/sec along 
the spiral trajectory of 15% overlap to sweep over the whole 
uncertainty cone for 55 sec. 

From Figs. 11 and 12, one can see that the scan time for 
spiral trajectory can be better than the other two trajectories. 

As pointed out previously, the square trajectories spend much 
more its scan beyond the uncertainty region. Similarly, 
the hexagonal trajectories may also spend some of its scan 
beyond the uncertainty region in order to ensure a complete 
coverage of the region. In general the path of trajectories 
starting from the center of uncertainty zone to its edge is 
shorter for spiral types .han for square and hexagonal types. 
Therefore, the spiral trajectories are recommended for Shuttle 
antenna scanning. 

Note that the constant speed along a trajectory yields a 
nonuniform angular velocity. In order to reveal the relation 
between the angular velocity and the constant speed along a 
trajectory, we plot the instantaneous angular velocity for a 
constant speed along a spiral trajectory, say 5°/sec, in Fig. 13 . 
Clearly it shows that in the neighborhood of the center of the 
ucnertainty cone, the Shuttle antenna is required to rotate 
more than 4 cycles/sec, which may appear too high for the Ku-band 
receiver antenna. To prevent this situation, the combination 
of the constant angular velocity around the center of the 
uncertainty zone and the constant speed along a trajectory for 
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outside area is proposed as the third alternative for the study. 

For the third alternative, an additional parameter, called 
switch point 0^, defining the point to change the motion of the 
Shuttle antenna from the constant angular speed to the constant 
speed along a trajectory, should be determined by the consideration 
of optimal dwell time for search cells. Here 'optimal' dwell time 
means to provide maximal available dwell time for signal detection. 

Thus, we shall discuss this problem in the next section. 

(3) Dwell Time - Mechanical and Electrical 

There are two kinds of dwell time involved in the acquisition 
system of spatial search— mechanical and electrical dwell times. 

The mechanical dwell time fs defined to be the time interval for 

Shuttle antenna to sweep along a trajectory from one search cell 

to an adjacent cell; while the electrical dwell time is the time 

interval for a signal energy detector to integrate the received signal energy 

and making a decision. Surely the electrical dwell time governed by 

the mechanical dwell time. In this study, we assume that the 

electrical dwell time is much smaller than the mechanical dwell 

time and that the electrical dwell time will be set uniformly for 

all search cells. The relation between these two dwell times is 

given as follows: 

where 


Tjji^ = mechanical dwell time 
Tjj = electrical dwell time. 
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For a constant speed along trajectories, the mechanical dwell 
time is uniform. Eq. (16) becomes 
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For a constant angular velocity, the mechanical dwell time varies 
depending on ti j location of cells. When the trajectory is near 
the center of the uncertainty zone, the number of cells swept by 
the Shuttle antenna per cycle is less and the mechatncal dwell 
time is longer, iherefore, the shortest mechanical dwell time 
available occurs at t'le outward area of the uncertainty cone. One 
may estimate the shortest mechanical dwell time by 






TT0 


(25) 


M 


where 

0 |^ = radius of uncertainty cone 
200 = beamwidth of antenna 

U 

f = revolution/sec 
A/ 20 g = overlap 

Next we discuss the case where the Shuttle antenna moves with constant angular 
velocity around the center of the uncertainty area then changes to 

a constant speed along a trajectory at a switching point 0^. We 
have plotted curves, as shown in Fig. 14, for kp = 0.1, and fixed 
total scan time = 55 sec. When the switching point at 0 .= 6°, 

the intersection of curves A and B yields the optimal electrical 
dwell time available, = 22.5 ms if one chose f^ = 0.26 cycle/sec 
for rotational speed within the switching point 0 ^. The curve A 
represents the minimum dwell time for cells within the constant 
rotation zone, while curve B is the dwell time of cells outside 
the constant rotation zone. Based on these sets of curves, a 
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system designer may determine the optimal electrical dwell time 
and also proper rotation speed. The constant speed for the outside 
switch point can be set to be the same speed as that when the 
Shuttle antenna reaches t‘ie switch point, that is 
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III. AVERAGE SCAN TIME (MECHANICAL) 

Scan time is one of the important parameters in the design 
of antenna acquisition systems. It is functions of the 
mechanical characteristtics of the Shuttle antenna, the scheme 
of signal energy detection and the statistical position of 
TORS relative to Shuttle users. The analysis of the scan 
time is rather complicated. In this study, we shall resort to 
computer simulation to obtain the scan time. However, we shall 
discuss the scan time of the Shuttle antenna without taking into 
account the effects of false-alarm and miss detection. This 
means the scan time will be determined for a chosen antenna 
trajectory, type of motion of gimbol systems and truncated 
Gaussian distribution for the TORS position. 

To derive the average scan time for the Shuttle antenna to 
point to the TORS position, we shall use Gauss quadrature rules, 
for which a general expression of scan time T conditioned on a 
known TORS positions is required. 

For the three proposed scanning trajectories the scan time 
T will be expressed in terms of either a rectangular coordinate 
system for spiral trajectories or in a rectangular coordinate 
system for two other trajectories--square and hexagonal. 

Denote the TORS position to be at a ^®EL’“AZ^ , where "el 
is in elevation component (radical direction in the projected 

K 

plane) and in azimuth component (angular direction). Let 

x(t) and y(t) be the rectangular coordinates of the trajectory 
of the Shuttle antenna motion on a projected plane at time 
instant t. One can relate it to the polar coordinates (0^,02) 




i 1 


quality 
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as follows 


®i(t) = + y^(t) in Elevation Direction 

^ (27) 

a,(t) = tan“^ in Azimuth Direction (28) 


In the analysis, it is assumed that, to start the spatial 
search of the TDRS position, the Shuttle antenna points tc the 
designated center of the uncertainty region of TDRS, say 
o(t) = (0,0) as shown in Fig. 15. Then the time required for 
the Shuttle antenna to scan up to the TDRS position, along 
the assigned scanning trajectory can be expressed as follows: 


T(a) = 


= f 


)ds along trajectory P. 

V 


where P. is the assigned trajectory of the Shuttle antenna motion, 

V 

He 

as seen in Fig. 15. The variable {^) is the velocity of the 
antenna motion along the trajectory at any instant t. More often, 
it is rather convenient to decompose the speed of antenna motion 
into two orthogonal components, either in polar form or in 
rectangular form. From geometric considerations one can simply 
express (^) in tne following forms: 

<E> ‘ J(“ (50) 


(~)^ + in Rectangular Form (31) 


Therefore the conditional scanning time T(a) on the known TDRS 
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TRAvTECTORY PKOJECTEn ONTO A PLAK'E PERPENDICULAR TO THE 
DIRECTION OF DESIGNATED CENTER OF UNCERTAINTY ZONE 
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position become 


T(a) = 




along Pj 


( /S\- 

(“l dt/ ^ \ dt/ 


for Polar Form 
(32a) 


T(5) = 


riU) 

ho,o) 

along Pj 




for Rectangular Form 


(32b) 


One can evaluate the scan time using the above formulation 
If the trajectory Is fixed and the TORS position Is known. 
Unfortunately, In practice the TORS position Is only known within 
a certain region. To characterize this unknovm, a truncated 
Gaussian distribution function p (a^ ,02) Is assumed for the TORS 
position. That means a Is a random variable with probability 
distribution p {0^,02)- With this uncertainty of the TORS 
position, the average scan time can be evaluated through 
Gauss quadrature formulas. By definition, the average scan 
time Is 


T(a) = E.[T(u)] 


= JJ' T(o^ 02)p(a^ .a2)do-| doj 


Here we shall assume that the variables and 02 are statistically 
Independent of each other, as Is the case In practice, I.e. 
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P( 0*1*02^ ~ P|(ft^)P2(®2^ 


Then one has 


. -HI 2 It 

T(a) = P/“l) f T(a^ * 012 ) 1 ^ 012)^32 

^ A JfS 


where 


0 ^ = the radius of the uncertainty region of TORS position 


The above double integral can be evaluated approximately 
through Gauss quadrature formulas as below: 


T = 


N, 

E E “U"2j<“U’“2k> 

4=1 k=l 


( 34 ) 


The sets of parameters and (‘02|(.“2k^ obtained 

from moments of random variables. A set of these parameters 
for the assumed probability distribution function (as shown in 
previous sections that the random variable is truncated Rician 
distribution and uniform distribution) is given in Table 1 
for = N ’2 = 5 with required moments of random variables. 

Through this scheme, the average scan time with two different 
types of the antenna motion along spiral trajectory of 30« overlap is 
computed vs o^, the standard deviation of the Gaussian distribution function 
The Shuttle antenna will take 31.5 sec to sweep over the whole 
uncertainty area if it is rotated with a constant angular 
velocity of 0.3 cycle/sec. The average scan time for this 
constant angular velocity is plotted in Fig. 16. For a mixed 
type of motion-'Sweeping with a constant angular velocity f^ 
then constant speed along the trajectory \>^ at a switching angle 0 


41 


'oClnG^m 


i 


! 4 


\ : 


j ^ 


I 


- V 


. >. ^ f 


jCinCc 


om 


Table 1. Moments and Parameters (w. ,a. ) of Gauss Quadrature 


Formulas. 


ir u' 


(1) Truncated Rician Variable (e^ = 1.0®, Op = 0.25®) 


(I-l)TH Central Moments 

1 l.OOOOOOOE+00 

2 3.1307825E-01 (Mean) 

3 2.6646434E-02 

4 2. 6520391 E-03 

5 2.2308420E-03 

6 6.2739134E-04 

7 3.5791496E-04 

8 1.5581140E-04 

9 8.4715485E-05 

10 4.4132538E-05 

11 2.4764277E-05 

(2) Uniform Random Variable (e = 

m 

(I-l)TH Central Moment 

1 l.OOOOOOOE+00 

2 5.0000000E-01 (Mean) 

3 8.3333333E-02 

4 0. 

5 1.2500000E-02 

6 0 . 

7 2.2321429E-03 

8 0 . 

9 4.3402778E-04 

10 0 . 

11 8.8778409E-05 


Gauss Quadrature Formula 


0 -2.4194992E-01 1 .2466771 E-01 

1 -9.5548526E-02 4.1508332E-01 

2 9.9011891E-02 3.6030230E-01 

3 3.2591406E-01 9.3469248E-02 


4 5.6917592E-01 6.4774140E-03 


Gauss Quadrature Formula 


0 -4.5308993E-01 1 .1846344E-01 

1 -2.6923466E-01 2.3931434E-01 

2 -7.9936058E-14 2.844444E-01 

3 2.6923466E-01 2. 3931434E-01 

4 4.5308993E-01 1 . 1846344E-01 
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the Shuttle antenna will take 61.6 sec to complete the scan over 

the uncertainty area if = 0.3 cycle/sec, = 5 044 degrees/sec. 

Its average scan time is also plotted in Fig. 17. Kjte that the 

unit of is the same as that in elevation unit/sec. 
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I V . KU-BAND RECEIVER— SIGNAL ENERGY DETECTOR 
(1 ) Receiver Model 

In this study tne receiver for spatial search is assumed 
to be a noncoherent signal energy detector, as shown in Fig. 19. 

The signal energy f>^om a square-law detector is used for detect- 
ing whether the Shuttle antenna is pointing toward the TDRS. 

During the spatial search, the Shuttle antenna is sweeping 
across the uncertainty region of the TORS along the assigned 
trajectory. At the same time the Shuttle receives signals 
transmitted from the TDRS. The signal power received by the 
Shuttle depends upon the offset angle between the boresight 
axis of the Shuttle antenna and the direction of the TDRS in 
the view of the Shuttle receiver. By comparing the received 
signal power with a threshold, one can determine that the 
boresight axis of the Shuttle antenna is pointing to the TDRS, 
therefore, stops the scanning of the Shuttle antenna. Denote 

the received signal to be S(t,e) plus the channel Gaussian 

2 

noise with zero mean and variance o^. "^he received signal may 
be expressed as follows: 

r(t,0) - S(c,e) + n(t) (35) 

where 

S(t, 0 ) = G(o) Sq(t)d(t) cos $(t) 

P = Maximum received signal pov/er when the boresight of the 
Shuttle antenna is pointing at the TDRS position 

Sq(t) = PN code waveform 

d(t) = Data waveform 

4>{t) = uipt + 
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Carrier frequency 
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i|»Q = Initial phase angle 

G(e) = Loss function of received signal due to offset 

of the boresight of the Shuttle antenna from the TORS 
position as defined in eq. (8) 

n(t) = A Gaussian random process with zero mean and 
variance 

n 

= n cos $(t) - n sin 4>(t) 

"c 1 

> = Baseband Gaussian random, processes with zero 
n I ^ 

^s means and variances 

n 


The variance of the channel noise is determined by the 
bandwidth of the IF filter Bjp and system noise temperature 
(noise spectral density Nq). Since the received signal i: 
wide spread in frequency, the IF filter must be designed to be 
wide enough to pass the signal. Therefore, the received signal 
power consists of not only the transmitted signal energy but 
also channel noise energy. When the available signal power 
at the Shuttle receiver is weak, the decision making based on 
the output of energy detection will be effected by the channel 
noise power; however, one can set proper threshold against the 
channel noise so that the probability of detecting = 0 
is reasonably high. Yet, when the available signal power is 
strong, the criterion of setting threshold level against the 
channel noise will cause false acquisition due to sidelobe 
of the Shuttle antenna, since the sidelobe of the Shuttle antenna 
is only 17 dB weaker than the strength of the main lobe of the 
antenna. In order to prevent this sidelobe acquisition, special 
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schemes must be implemented. Some recommended techniques for 
avoidance of sidelobe acquisition will be discussed separately. 

The signal detection for the Ku-band receiver can be 
formulated as follows: 

Denote the output of ine square law detector to be 

y(t.e) = r^{t,e) 

depending on the offset angle between the boresight axis of the 
Shuttle antenna and the direction of the TDRS. The output of the 
integrator is 


(36) 


Z(6) 


■rj 


r'^(t,e)dt 


(37) 


The detector is to determine one of two hypotheses: 


Hq*. Z(e) 
Z{0) 


|el 16a 
|ei < 


(38) 


where 6^= allowable offset angle of the Shuttle antenna. The integration 
of Z(e) will be approximated by the sum of N samples taken at 
a Nyquist rate T^, which is equal to Bjp. 


Z(e) = i 2] r^(kT^.e) 
n=1 

where M = integer part of BjpT^. The statistical characteristics 
of the variable Z(e) can be determined from that of a single 
sample r(kTg.e), to be discussed in the next section. Here we 
shall further make cne assumption that e. ^ 0. So the hypotheses 

a 


M 


(39) 


should be modified as 
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Hq: Z(e) 
H,: Z(e) 


for 0 > 0 


for 0 = 0 


The following analysis can be modified for small values of 

O 


(2) Statistical Characteristics of Decision Variable 

By the assumption that n^ and n^ are two lowpass Gaussian 

2 

random variables of zero means and variance o^, the probability 
density function of tne envelope R and phase ^ of the signal r(t,9) 


is well known 


/2PG^(e) + - 2v7PG(0)S„(t)d(t)Rcos ^ 


P(R.<.l0,Sq.d) =< " 

0 for R < 0 


where 


R = [v^G(e)S (t)d(t)+n^+ n^ 


for R > 0 


^ = tan 
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Define 


i^(e) = = input signal-to-noise ratio (SNR) 


P(R.i{»|0,S^.d) = 




Sq(t)d(t) cos \l> 


for R > 0 


R < 0 


Averaging P(R,i{)!0) with respect to one has 


$ 
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P{Rlo.S^.d) = 


Assuming that 


exp 




R < 0 


Sq(t) 


with equal probabilities 


R > 0 


1 

d(t) = 1 } with equal probabilities 


then one has 


f*(rJo) 


where 


(- d„) - »'('')) 'o )*'<’>) 


for r^ < 0 


r = sample of R at t = iT 
* ^ 
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Now we can derive the probability density function of Z(o) from 
P/r^jo) through tiie characteristic function approach, which can 
be found in [1]. Omitting the detailed derivation, we express 


the f'esuU as 


P(Z|0) = 


(“■'<"> Xfp) •'O'- ^ If 

n t , 


I < 0 
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When e = 0, the probability of detecting the presence of 


TORS can be computed 

Pj(0=O) = 


where Th is a threshold level of the detection. Defining - 
2 

MZ/2o^, one has 

Pd(e=0) = q„(hp^(0),^) (47) 

\ / 

where 

Q|^(x,y) = J (|)”“ exp(-z-x)In^ .j( 2 /)^)d 2 

y 

The function Qj^(x,y) is known as generalized Marcum's Q-function. 

When the offset angle o is nonzero 0g), the probability 
of false alarm can be expressed as 

ORIGINAL PAGE IS 
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= QMfMP^(o), ^ (48) 


as 


J P{Z|e=0)dZ 


Th 


,2o^P^(0) 


■C 4 ( 

-•{"■'"S') 


\M-l/2 / 

) exp 


MZ 


2a: 


dZ 


(46) 







om 

which is a function of the offset angle 0. The criterion for 
choosing Th is the prevention of false detection caused by 
the second sidelobc of tne Shuttle antenna. In practice, for 
a given false alarm probabil ity.one needs to establish the 
threshold level Th. However, solving for Th from eq. ( 48) 
for a given value of appears to be nontrivial unless s^-e 

approximation is made. Thus one may assume that p (oj'O) 0. 


Pf,(0i<0) 


1i»i Q„(Mp^(e), -5 ) 


P (0j'0) >0 

MTh 
M-1 


exp(-2)dZ 


/MTh\"'^ / MTh\ ty 

(>?) 


1 


n' k=0 


(M-k 




M 7 h\ir" 

- 2‘ 

n 


for M > 1 (49) 

Using iterative scheme one can solve for Th for a given value 
P^^(o/0). In fact the above expression can be rewritten as 


Pf(pj'O) =1-1 


Ky ■ ”■’) 


(50) 


where 


I(u,s) = incomplete gamma function 


/ 


S [ . 

«_e><£(,vx.) dx 


S! 


The function I(u,s) has been tabulated 


Thus one may 


*Pearson, K. , Ta^ijcs of jncomplcte r Cambridge University 

Press, 1934. 
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conveniently solve for the threshold Th. Note that the approxi- 
nation p (o^O) ^ 0 is valid when the noise variance is of the 
order of the first sidelobe magnitude (about SNR < 17 dB). Other- 
wise, alternative approximation may be taken. When signal-to- 
noise ratio SNR is high and the value of M is large, we shall use 
an equivalent Gaussian distribution for P(Zlo). 

, / (Z-m,(6))^\ 

° J77 "P r ■■ 2oj,(e) “ j 

whdre 

012 ( 0 ) = M(1 + p^(o)) (52a) 

0 ^( 0 ) = M(1 + 2p^(6)) (52b) 


Then the probability of false alarm can be simply expressed as 


where 



(53) 


Thus the threshold Th becomes 

Th = *2o^/M (54) 

To evaluate Th, we substitute 0 by e^p, the offset angle at 
which the first sidelobe of the antenna is located. The function 
4 >"^(') is tne inverse function of erfc(-) function. 

( 4 ) Threshold Levels of S ignal E nergy Detect or 

The choice of threshold level in signal energy detection 
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depend:> upon the specification of false alarm probability and 
probability of detection and carrier-to-noise (CNR) ratio. One 
may start to choose initially a threshold Th^ to meet the probability 
of false alarm due to channel noise. When CNR becomes larger, 
as indicated previously, one has to verify the chance of being 
acquiring the TORS with the sidelobes of Shuttle antenna. If 
it is higher than one needs to adjust the threshold level 
higher so that the probability P^^^ of detecting the TORS with 
sidelobe of antenna is less than specs. 

However, in practice the CNR is not known at a Shuttle 
receiver, it is better for hardware implementation that a 
uniform threshold level for all operating range of CNR is set 
up. For a Ku-band communication link, the operating range of 
CNR can be as low as 60dB-fl^ and as liigh asSOdB-Hz. With this wide 
range of CNR, (60 dP Hz - 80 dB Hz) one \.ill find it impossible 
to sot a unifom threshold level without having false acquisition 
due to sidelobe of the Shuttle antenna. This can be clearly 
seen when the signal energy at the output of detector is plotted 
for the main lobe against the sidelobe of Shuttle antenna, as 
shewn in Fig. 19. Note that the signal energy received by a 
sidelobe of Shuttle antennas at CNR = 80 dB-Hz, can be 
stronger than the signal energy received by the main lobe of 
the antenna at CNR = 60 dP-Hz. Therefore, the threshold set 
up for CNR - 60 dP-Hz will be crossed by the signal energy of 
the sidelobe at CNR = 80 dP-Hz. 

In order to avoid the false acquisition due to the sidelobe 
of the Shuttle antenna, a second threshold should bo used to 
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gure 19. Received Signal Energy. 

•IF Bandwidth 6 MHz 
•Integration Time 6 
•Filter Loss 0.5 dB 
•Ohter Loss 1.0 dB 
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distinguish the main lobe acquisition from a fi>'St sidelobe 
acquisition. Here it is proposed to use the second threshold 
level at CNR = 80 dB-Hz to guarantee the specs of the probability 
of false alarm due to the sidelobe for high SNR. 

One also should note that for a wide operating range of 
CNR (60 dB-Hz to 80 dB-Hz), the possibility of having a second 
sidelobe acquisition occurs if the second sidelobe of antenna 
pattern is not lower than the main lobe by about the range 


of CNR. 




V. ACQUISITION SYSTEM OF KU-BAND RECEIVER 


1 . Three Acquisition Strategies Studied 

The acquisition system proposed for Ku-band receivers 
consists of two major parts— coarse spatial acquisition and 
autotrack >cquisition. Here, we shall confine ourselves to 
the study of coarse spatial acquisition system. From the 
system point of view, it is expected that at the end of coarse 
acquisition, the boresight axis of the Shuttle antenna must 
be at least aligned with the direction of the TORS position to 
within the pull-in range of the autotrack system. Three 
acquisition algorithms studied here can all acquire the TDRS 
within the main lobe of the Shuttle antenna when the received CNR 
at the receiver is not too high {< 75 dB-Hz). However, with a wide 
dynamic range of CNR, one can identify the acquisition ranges of 
these three schemes to be within sidelobes of the Shuttle antenna, 
within the nulls of the main lobe and near 3 dB beamwidth of the 
antenna, respectively. 

The acquisition schemes used for the coarse spatial search 
can be described as below; 

(1) Normal scan - The boresight axis of the Shuttle antenna 
is offset at most from the direction of the TDRS position 
by an angle equal to the range of the first sidelobes 
of antennas. In this case the hit declared by the 
acquisition system may be due to either the main lobe or 
the first sidelobe of the antenna and this is dependent 
on the signal -to-noise ratio. Thus the autotrack 
system should be initiated to search for the TDRS position 
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by the main lobe of the Shuttle antenna. 

(2) Normal scan plus the avoidance of the sidelobe acquisi tion- 

The scan will put the boresight axis 

of the Shuttle antenna within null angle cf the 

main-lobe of the antenna from the TORS position. 

For this case, the acquisition procedure is ended 

at the acquisition of the TORS position with the 

mainlobe of the antenna. Thus the range of the 

autotrack operation will be much smaller than the 

previous case. 

(3) Normal scan plus fine-search - At the end of the 
acquisition, the boresight axis of the Shuttle 
antenna will be within 3 dB beamwidth of the antenna 
from the TORS position. This scheme shall not only 
avoid the sidelobe acquisition but also give a finer 
range for the autotrack acquisition. Surely the 
acquisition strategy for this case is much more complicated than 
the previous cases. 

2. Acquisition Strategies 

The acquisition strategies for the three acquisition schemes 
are depicted by flow diagrams in Figs. 20 ,21 and 22 , respectively. 

For the normal scan, there are only two stages within the 
acquisition— routine scan and verification. While the Shuttle 
antenna sweeps through the whole uncertainty area along a 
designated trajectory, the received signal energy at the 
output of the signal detector is compared against a preset 
threshold THl . If the threshold is not exceeded, the search 
continues until the end of the whole area. If this is the 


case, the search fails to acquire and the search should restart.^ 


o/u 


I I 


I V 





Figure 20. Acquisition Strategy with Normal Scan. 
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On the other hand, if the threshold is exceeded, a hit is 
claimed and verification stage is entered. Once a second hit is 

obtained, the acquisition of the TORS is declared and the procedure 

of the spatial search s stopped; otherwire, the search continues. 

The second acquisition scheme is designed to exclude +he 
side!, 'e acquisition which occurs for large values of CNR. Two 
threshold levels are set up to distinguish the sidelobe acquisitions 
from the mainlobe ones. The following rules for two threshold 
levels are: 

(1) For a specified probability of false-a1arm, and 
an operating range of CNR, a threshold THl is 
computed at the lowest value of CNR. 

(2) At the largest value of CNR, the second threshold TH2 
against the detection due to the sidelobe r antennas 
is computed for the specified 

(3) The threshold level THl should be used initially in the 
stage of acquisition. Once a hit is indicated, threshold 
TH2 is used for a mini scan to verify any sidelobe 
acquisition. 

The status of the second acquisition can be briefly 
classified into three stages--routine scan, mini-scan ana 
verification. The procedure of the acquisition is similar to 
the normal scan, except that v.'hen a hit is indicated, the 
acquisition enters a mini -scan stage. The purpose of the niini- 
scan is to verify whether the hit H in the routine scan is caused 
by the main lobe of the aiitenna or by its sidelobe. The 
trajectory of the mini-scan is a circle centered at the position 
of the hit H with radius R = 2.5® (for Og = 0.8 degree of antenna). 
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which is the separation distance between the first sidelobe and 
the mainlobe, see Fig. 23. Regardless of hit or not after exhausting 
the mini-scan trajectory, the acquisition enters a verification 
stage, as indicated in Fig. 21. 

There are two different cases for verification. First, the 
verification of a hit during mini-scan searcn requires nc change 
in threshold level. If a second hit is observed during this 
verification stage, the search is terminated and the system 
enters autotrack mode for refine acquisition; otherwise, it 

t 

returns to mini -scan. For t?;e second case, in which no hit was 
observed in mini-scan, the antenna is moved back to the position 
of hit H in the routine scan for a verification. This case will 
lower the threshold level to THl. If a hit still is observed, 
the acquisition procedure terminates and enters the autotrack 
mode; otherwise, the routine scan is resumed. 

In addition to the strategy used in the second scheme, 
the third scheme will scan the cells adjacent to the hit H in 
routine scan to acquire the TORS position within the beamwidth 
if the hit was caused by small offset from the TORS position (it 
happens when the CNR is around 80 dB Hz). The signal energies 
received in the adjacent cells are compared and the cell with 
the strongest signal energy will be chosen as the right hit. 

If there is no hit in the adajcent cells, the system will 
proceed the mini -scan as proposed for the sect d scheme. The 
flow diagram of tivj thirc .. neme is shown in Fig. 22. 

( 3 ) Performance of Acquisition Systems 

The performance of the acquisition systems will be measured 
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by the average acquisition time in acquiring the TORS position. The 
mathematical analysis of the acquisition performance is 
rather complicated because of the following reasons: 

(1) Nonuniform probability distribution of the TORS 
position. 

(2) Nonunifom! mechanical dwell time for each cell. 

(3) Complicated acquisition strategies, such as mini- 
scan and fine scan. 

However, a computer simulation program was developed for the 
study of the performance of acquisitions for various system 
conditions such as the probability distribution of the TORS 
position, trajectories of Shuttle antennas, motion along the 
trajectory, CNR, probability of false alarm, probability of 
detection and different acquisition schemes. Figs. 24-26 
give the performance of typical acquisition systems obtained using 
the simulation program. The simulation program and its capabilities 
are documented in Volume II of this report. 

The scan time required to acquire the TORS position was 

obtained through Monte Carlo approach. There are 500 simulated 

TORS positions for each run. The analytic results given in 

Figs. 20 and 21 have been closely verified through the simulation 

program for CNR = 60 dB-Hz and normal scan acquisition algorithm. 

In addition, the third acquisition algorithm with combination of 

constant angular velocity and the constant speed along a spiral 

trajectory was also tested. Their results are plotted in Fig. 26 

(also in Table II) for CNR = 60 dB-Hz and CNR = 80 dB-Hz. In 

these cases, we assume that the transition from a normal scan 
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Figure 24. Average Spatial Acquisition Time for Constant 
Angular Velocity. 
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Figure 26. AVERAGE SPATIAL ACQUISITION TIME FOR THIRD 

ACQUISITIOri STRATEGY AND A MIXED MOTION ALONG 
SPIRAL TRAJECTORY OF 30S OVERLAP. 
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to adjacent scan, or to sidelobe scan and vice versa, takes 0.0/sec 
and that motions along both the adjacent scan and the sidelobe scan 
are the constant angular velocity with respect to their scan center. 
It shows that at high CNk, the spatial acquisition is much shorter 
than that at a lower CNR, since the sidelobe scan and adjacent scan 
help shorten the scan time. Therefore, the third acquisition 
algorithm should be recommended for the consideration of the Ku-band 
acquisition system if the operating range of the CNR is wide, 
such as that used for this study, and if the side-lobe acquisition 
imposes a potential problem to the system performance. 
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Table II. Acquisition Time for Third acquisition Algorithm 
with Mixed Motion Along a Spiral Trajectory of 
30%. 


#*e, 


m 



60 dB-Hz 


Average Scan Time 

1.5 

11.48 

2.5 

15.4 

3.5 

23.1 

4.5 

22.75 


= 10% 0g = 1.6 

= 5.5®, f^ =0.3 cycle/sec. 


80 dB-Hz 

Average Scan Time 
2.675 
5.386 
9.287 
12.58 


= 5.044 degree/sec 
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VI. COMPUTER SIMULATION OF KU-BAND ANTENNA POINTING 

ACQUISITION SYSTEM 

(1 ) Introduction 

The Ku-band antenna pointing acquisition systems based on 
the noncoherent signal energy detection and scanning scheme 
discussed previously have been simulated by digital computers 
for the tradeoff study in the design of the Shuttle antenna 
acquisition. Several acquisition algorithms implemented are 
given to provide flexibility in accommodating the different 
operating ranges of autotrack systems. The acquisition strategies 
are primarily oriented toward the avoidance of the first side 
lobe acquisition in the spatial search. In addition, the 
simulation also provides various options on the scan trajectories, 
scan rates and a wide range of CNR. It essentially integrates 
the effect of interaction of spatial search acquisition schemes 
and the noncoherent signal energy detection. So the simulation 
program is a^;le to support the design of hardware development. 

The simulation program was written in Fortran IV for the 
Univac 1100 series computer. Origina’ly, it was developed for 
time-shared operation, which is available at LINCOM Corporation. 
However, the program can also be used for batch processing with 
some modification. The program consists of 19 internal subroutine 
references and two external subroutines MERFIC (from IMSLS package) 
and BSSL (from MATHPACK) which are normally available for 
computer systems. 

The computing time to run the simulation program depends on 
the number of tests on the TORS position and the number of 
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samples taken for channel noises. These two variables are also 
provided as options for a user. The higher the values of these 
two variables, the higher the accuracy in predicting the acquisition 
performance and the longer the CPU time. The CPU time for a 
typical case is well within 10 sec if the number of noise samples 
is one and the number of simulated TORS position is less than 500. 

The simulation program provides various options for users 
to test system performance for a wide range of system parameters. 

The following options included are: 

(1) The size of uncertainty cone of the TORS position and 
its statistical distribution function over the cone. 

(2) The antenna pattern of the Shuttle receiver (3 dB 
beamwi dth ) . 

(3) The types of antenna scan trajectories--square, hexa- 
gonal or spiral with different overlap. 

(4) The scan rates— constant angular velocity, constant 
speed along a trajectory or combination of both previous 
cases (also a switching angle 9^). 

(5) The probability of false-alarm and the operating range 
of carrier-to-noise ratio (CNR). 

(6) The receiver parameters--dwell time of the signal energy 
detector, IF filter bandwidth, insertion loss and other 


circuit losses. 

(7) The number of TORS positions used in the Monte Carlo 
simulation. 


(0) The number of samples taken for channel noises. 
(9) The types of acquisition algorithms. 
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(10) The choices of a detailed printout or final results 
in the simulation. 

The simulation program was structured as functional modules 
so that each subroutine can be easily replaced to adopt to 
any required function. For instance, the antenna pattern used 
in the simulation can be simply replaced by other types of 
antenna pattern. To do so, one should remove the subrou-.ne 
ANTNNA and substitute an equivalent subroutine ANTNNA. 

The simulation results are presented in two versions--one 
for a time-shared mode and the other for a batch process mode. 

For the time-shared mode, the computer asks for input data 
through interaction with users. A typical output generated by the 
simulation program consists of the following important parameters: 

(1) The threshold THl and TH2 used for the signal energy 
detector in a normal scan and in mini scans (such as 
adjacent and sidelobe scans), respectively. 

(2) Total number of cells to cover the whole uncertainty 
region. 

(3) Maximum scanning time (= time spent to reach the last 
cell through normal path only). 

(4) Average scanning time and corresponding standard 
deviation. 

(5) "Probabilities" of detection, miss and false alarm. 

A sample printout is given in Fig. 27. 

(2) Functional Diagram of S-»’mulation Program 

To clarify the main functions performed in the simulation, 
a flow diagram of the simulation program is depicted in Fig. 28. 
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RADIUS OF UNCERTAINTY REGION: 10.0000 DEGREES 
BEAKWIDTH; 1.6000 DEGREES 
SPIRAL TRAJECTORY: OVERLAP = IC.OOOO ? 


CONSTANT ROTATIONAL SPEED: 
PROBABILITY OF FALSE ALARM: 
CNR EXPECTED RANGE: 
INTEGRATION TIME: 

IF FILTER BANDWIDTH: 
INSERTION LOSS: 

OTHER LOSSES: 


108.0000 DEGREES/S 
. 1000-02 

60.0000 DB - 80.0000 DB 

5.0000 MS 

6.0000 MHZ 
1.5000 DB 

.0000 DB 


THE PROBABILITY OF FALSE ALARM WAS RAISED TO 0.1 TO AVOID SIDELOB 

. . .E DETECTION 
WHEN CNR = 80.00 DB/HZ 


THRESHOLDS: TH1 = 1.0568141 TH2 = 1.2276665 

NUMBER OF CELLS: ?02 

MAXIMUM SCANNING TIME = .3217689+02 SECONDS 
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NUMBER OF SATELLITE POSITIONS GENERATED: 500 

STANDARD DEVIATION: .1000000+01 


CARRIER-TO-NOISE RATIO: 

ACQUISITION AVERAGE TIME: 

STANDARD DEVIATION: 

PROBABILITY OF DETECTION: 

PROBABILITY OF DETECTION (SL. VAR.): 

PROBABILITY OF MISS: 

PROBABILITY OF FALSE ALARM (IN THE SIMULATION): 


.6000000+02 DB 
.3023562+01 S 
.2362385+01 S 
. 1000000+01 
.0000000 
.0000000 
.0000000 


Figure 27a. A Sample of Input and Output Data. 
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Figure 28. Functional Flow Diagram of the Simulation Program. 
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The TORS position is simulated by a random generator with a 
specified distribution function. To measure the spatial 
acquisition time of the TORS position, a set of simulated TORS 
positions is tested for a specified trajectory and scan rates 
(speed for constant angular velocity and/or constant velocity 
along a trajectory). The simulated data reported here are all 
obtained by 500 simulated TORS positions. 

In the simulation program, the acquisition monitor is not 
a separate unit, but -mbedded in the subroutine ACQ. This is 
inq)lemented just for the convenience in programming. The 
detailed structure of the program is referred to Volume V. 

(3) Capabilities o f ~^ ftware Package 

The program devtis-ped for the Ku-band antenna pointing sy'-tem 
provides various options in the ;ize of uncertainty cone for the 
TORS positions, the type of scanning and trajectories, and spatial 
acquisition algorithms. Here are the capabilities of the software 
p£ ckaqe : 

(1) Conical angle of the TORS position can be 4®, 8' or 10®. 

2 

(2) The variance Op of a priori probability distribution 
function of the TORS position can vary from 0.1® to 10°. 

(3) Antenna pattern can have beamwidth > 10®. 

(4) Three scan trajectories and wide range of scan rates 
(but limited by electrical dv/ell time) are available. 

(5) The noncoherent signal detector can be operated over 

a wide rt.nge of CNR and BjpT product greater than 1000. 

(6) Spatial acquisitiun algorithms with/without avoidance 
of sidelobe acquisition are selectable. 
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(?) ProhaMlity of fvUsc alann ran be from O.f' down to 10*^. 

(B) The rumiber of TORS positions in Kmte Carlo si iulalion 
can be selected by users. Rut it is limited by the CPU 
time in executing the simulation program. 

In addition, the software package also provides a useful 
feature for system designers to interact with a co.:;puter if the 
time-shared facility is available. Upon user's choices, a detailed 
spatial acquisition procedure ran be printed out for a closed 
examination on the behavior of acquisition systems, finally 
one should also note that the application of the software 
package should not be limited to the Ku-band system I'aranscters 
even though it is oriented tow.u'd the study of Ku anter.na pointing 
systems. 

(4 ) Coni|niter Preyram Ut i 11 rat ion 

The simulation progi’am was ..csigned for two nodes of usago-- 
batch process and tine-shared modes. Since the computer proipwn 
was originally developed through time-shared node, it provides the 
mechanisn of interactions betwet-n computers and users. However, this 
set of interactions should be suppressed for the batch processing. 

In order to run the program in batch, one should provide the samt' set 
of input data as one type in input data in the tine-shared mode, 
except for the first input data. To illustrate the details of the 
input data required, several examples are given in figs. CV and CT 
to show the options available in the simulation and the cof'rospond’og 
sets of input data for tf.o batch process are also listed. 
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(5) Func tion al Descr ipt ions of Subr outines 

The functions of each subroutine used in the computer program 
for UNIVAC 1100 series are briefly described one-by-one here. 

Their flow charts and detailed documentations are given in 
Volume II separately. The description of subroutines is 
given in alphabetical order to be consistent with the computer 
printouts from UNIVAC llOS time-shared mode. The program is, 
in fact started from MAIN then SIMULA. The subroutine SIMULA, 
the core of the simulation, cells and monitors the operations 
of the Monte Carlo simulation. 

Subroutine ACQ 

This simulates the acquisition algorithms discussed 
previously. The variable NULL, used to select the type of 
the acquisition algorithms for the simulation, is defined as 


follows: 


NULL = 0 for normal scan 


1 for normal scan plus sidelobe scan 

2 for normal scan, adjacent scan and 
sidelobe scan 

In the process or the spatial acquisition, the time spent in 
any scan, including the transition from one scan to another 
and verification time, is accumulated as one of the outputs of 
the subroutine. The other important output parameters of the 
subroutine are: 

(1) The status of the spatial soarch--success or failure 
(IPMISS and IPDET). 

(2) How to acquiro--the detection was declared by a normal 
scan orplus sidelobe scans (IPDSL). 
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(3) The acquired location of the TORS (KHIT (IHIT)). 
Subroutine ANTNNA 

The antenna pattern of the simnlated Shuttle receivers is 
assumed to have a first sidelobe of -175 dB from the peak of the 
main lobe. The antenna gain is computed from the offset angle 
between the simulated TORS position and the boresight axis of 
the Shuttle. Hence, the offset angle (TH) is an input variable 
in addition to the parameters characterizing the antenna (THB, 
PATTC). A control variable K is used for computing the 
antenna parameters as shown here 


K = -2 

for computing the angle of first null 

-1 

for computing the angle of the peak of 
of the second sidelobe 

0 

for computing the angle of the peak of the 
first sidelobe 

+1 

for -3 dB beamwidth 

2 

for normal ge in computing the gain 

of the antenna at offset TH. 

Subroutine CELL 



For a specified trajectory (KSCAN), the subroutine computes 


the locations of search cells (XCTR.YCTR) along the trajectory 
and the scan time to reach each search cell from the designated 
center of the uncertainty zone of the TORS position. The scan 
time is computed by calling the subroutine TCENTR. Hence, it 
requires the following input variables 

KSCAN = 1 for square trajectory 

2 for hexagonal trajectory 
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MOTION = 1 for constant rotational speed 

2 for constant speed along trajectory 

3 for constant rotational speed first 
then constant speed alc;;g trajectory 

VANG = angular velocity for constant rotational 
speed (in radian/sec) 

VLIN = constant speed along trajectory (in radian/sec) 
To ensure the electrical dwell time to be smaller than times 
of the minimal mechanical dwell time, a warning is printed out 
if the antenna scan is too fast. In the subroutine is set 
to be 1J3 . 

Subroutine GAUSS 

This is a Gauss random generator whose mean and standard 
deviation are specified by AVR and SIGMA, respectively. The 
Gaussian samples are obtained in pair by the following simple 


relations 


= o(-2.1n x^) cos 2 uX 2 + v 


1 /2 

y^ = o(-2.1n x^ ) ' sin 2 nX 2 + u 


where and X 2 are a pair of independent random variables uniformly 
distributed between 0 and 1, and p and a are its mean and variance, 
respectively. The algorithm has been tested for its mean, van’ance 
and skew coefficient. The results show that the alaorithm gives 
V high confide.’ in its statistical natu*'e. 

The signal energy detector is simulated by this subroutine. 

The received signal energy (sum of signal power and noise power) 
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is comparGd with a preset threshold (TH). If the threshold is 
crossed, a hit is declared (HIT=.TRUE. ); otherwise, no hit is given. 
The input variables are the position of the TORS and the location 
of the boresight axis of the Shuttle antenna. The offset angle 
lETA is computed, then the antenna gain by calling subroutine 
ANTNNA, the signal power and noise power. 

Subroutine HPOSTN 

This subroutine computes the ceiiters of hexagonal cells along 
the hexagonal trajectory. The computation is proceeded from 
K— cell to (K+D-- cell. 

Subroutine ITR 

This subroutine implements the iterative algorithm needed 
in solving for eqs. (20) and (21). The control variable K 
is used to indicate the application of the algorithm to eq. (20) 
or eq. (21). 


The function of the MAIN defines the dimension of variables 
used in the simulation. The program has been set to deal with 
the maximal number of search cells to ^^e 500. If the rnimber 
of the search cells goes beyond 500, one should expand dimensions 
of all variables in the first dimension statement. 

MAIN/MAP 

This is a set-up for executing the simulation program. Two 
system subroutine packaaes--MATHPACK and IMSL--are called to 


satisfy the external references. 
Subroutine NOISE 
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The simulated channel noise power at the output of the integrator 
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is computed here. The number of samples for equivalent Gaussian 
noises is set by the variable KNOISE. 

Subrou ti ne QPOSTN 


It computes the center of (K+l)~ search cell from the known 
center of the K— search cell along the square trajectory. 
Subroutine SIMUL A 

This is the core of the simulation program. It reads in the 
data, prints the outputs and also monitors the operation of the 
program. The detailed operation of the subroutine is referred to 

^‘n Volume V. 

Subroutine SLCELL 

The search cells alons a sidelobe trajectory and adjacent 
cells to a hit cell are computed here. The adjacent cells are 
structured as those for hexagonal cells. Hence, the number of 
the adjacent cells to be scanned is fixed as six. However, 
the number of cells for the sidelobe trajectory varies depending 
upon the type of trajectory used for normal scan. 

Subroutine SPOSTN 


s t 

It computes the center of (K+1) — cell along a spiral 
trajectory from the known center of the K— cell. 

Subroutine TCENTR 

It computes the scan time required for the Shuttle antenna 
to sweep along a trajectory to reach a spatial point in the 
uncertainty zone of the TORS position. It requires to input 
t!ie type of motion, scan rate (constant rotation cnd/or constant 
speed along a trajectory) and the coordinates of the point. 
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Sub ro utine TORS 

The subroutine venerates a simulated TDRS position according 
to a specified Gaussian distribution. If a Gaussian sample falls 
outside the uncertainty zone of the TORS position, it is discarded 
and another sample is taken. Its input variables are SIGMA, 

THM2, and its output variables XTDRS and YTDRS. 

Subroutine TH1TH2 

The setup of threshold levels used for spatial acquisition 
is done in this subroutine. The threshold THl is first computed 
for the given probability of false-alarm under the channel noise 
condition. Then it checks the probability of detection 
due to the second sidelobe of the Shuttle antenna at the 
maximal value of CNR in the specified operating range. If 
is larger than the given probability of false alarm, the threshold 
is adjusted to yield at least 0.1 probability of false alarm 
due to second sidelobe at the strongest CNR. This implies that 
the probability of false-alarm at lowest CNR is much smaller 
than the specified one. 

The threshold TH2 is computed based on the assurance of 
having the speci -'ed probability of false alarm due to the first 
sidelobe at strongest CNR. 

Subroutine UNIFOR 

This is a uniform random generator over (0,1). A pair 
of output samples are generated for each cell. 
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The Ku-band antenna pointing system has been discussed and 
simulated by digital computers. Through the analytical and 
experimental results presented in this report, one may lead to 
the conclusion that the spiral trajectories with various scan 
rates (mixed motion) should provide the flexibility for the Ku- 
ba.id antenna spatial scan over the TORS uncertainty region and that 
the tnird acquisition algorithm proposed should not only avoid 
the side-lobe acquisition problem, but also yields a shorter 
acquisition time especially wlien the received signal at the 
Shuttle is strong. 

With the various optio"" ,-.'cii<ble in the software package 
developed for this study, -.xu. ' ould be able to perform, with ease, 
tradeoff study on the system pararnetet's to yield an "optimal" 
system per romance. 
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